Queuing transitions in the asymmetric simple exclusion process 
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Stochastic driven flow along a channel can be modeled by the asymmetric simple exclusion 
process. We confirm numerically the presence of a dynamic queuing phase transition at a nonzero 
obstruction strength, and establish its scaling properties. Below the transition, the traflic jam is 
macroscopic in the sense that the length of the queue scales linearly with system size. Above the 
transition, only a power-law shaped queue remains. Its density profile scales as 5p ~ x~" with 
= |, and X is the distance from the obstacle. We construct a heuristic argument, indicating 
that the exponent f = i is universal and independent of the dynamic exponent of the underlying 
dynamic process. Fast bonds create only power-law shaped depletion queues, and with an exponent 
that could be equal to v = ^, but the numerical results yield consistently somewhat smaller values 
V ~ 0.63(3). The implications of these results to faceting of growing interfaces and localization of 
directed polymers in random media, both in the presence of a columnar defect are pointed out as 
well. 
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I. INTRODUCTION 

Queuing is a common nonequilibrium phenomenon in 
nature. It appears, e.g., in stochastic-type driven trans- 
port through narrow channels, where applications can 
range from electron transport along nanowires to traf- 
fic flow on highways. Queuing has received a lot of at- 
tention from the theoretical side for over a decade [1,2]. 
It is well established that many driven flow processes 
belong to the same universality class as Kardar-Parisi- 
Zhang-(KPZ) type growth of one-dimensional (ID) in- 
terfaces [3,4]. In particular, the so-called asymmetric 
simple exclusion process (ASEP) [5] maps exactly onto 
the so-called body-centered solid-on-solid lattice version 
of KPZ growth [6]. This model has been used to de- 
scribe biopolymerization [7], gel electronics [8], directed 
polymers in random media [9], traffic jams [10,11], and 
the fluctuations of shock fronts [12-15]. In the case of 
periodic boundary conditions, the time development of 
the ASEP is exactly soluble by the Bethe ansatz [5]. 
For a wider class of setups, the exact stationary state 
has been constructed as well using the so-called matrix 
method [2,13,15-17]. 

The starting point and motivation for the study pre- 
sented here was actually not queuing in driven flow but 
faceting in KPZ growth. Slow flameless combustion of 
paper produces ID burning fronts that evolve in time 
according to KPZ- type growth. The results of early ex- 
periments [18] seemed to deviate from KPZ behavior, but 
more recent investigations demonstrated that for length 
scales larger than about 5-10 mm depending on paper 
structure, random pinning effects do not play a role 
any more, and that the interface obeys ID KPZ scal- 
ing [19]. In these experiments, the paper is impregnated 
with KNO3 to provide the oxygen source necessary for 



maintaining the slow-combustion process. The burning 
speed can also be controlled by the KNO3 concentra- 
tion. In particular, the rate can be enhanced or reduced 
in a narrow strip along the burning direction. The lat- 
ter experiments nicely illustrate the presence of nonlinear 
terms in the equation of motion; the burning front facets 
for enhanced concentrations but not for reduced ones. 
The details of these experiments, as well as the match- 
ing of the experimental data with our numerical results 
for the slow- and fast-bond ASEP will be published sep- 
arately [20]. Here we only present the ASEP queuing 
perspective. 

One of the fundamental questions in driven flow is 
whether a static obstruction, such as a slow bond, always 
results in a traffic jam, or whether stochastic fluctuations 
destroy the queue at weak obstacle strengths. Such a 
vanishing of the queue as a function of the slow-bond 
strength represents a dynamic phase transition. The size 
of the queue is the order parameter, i.e., it being flnite or 
infinite in length; or more precisely, whether the number 
of "cars" in the queue scales and diverges with the sys- 
tem size or remains finite (in obvious analogy with macro- 
scopic occupation of the ground state in equilibrium Bose 
condensation). The existence of such a transition, its 
scaling properties, the shape of the density profile near 
the obstruction, and also whether information percolates 
through the slow bond, are the most important issues. 

The ASEP is one of the simplest nonequilibrium driven 
dynamic processes displaying queuing phenomena. Parti- 
cles move stochastically along a chain of sites 1 < a: < iV^, 
only in one direction, with hopping probability p under 
the constraint that the particles can neither pass each 
other nor occupy the same site, = 0, 1. We use ran- 
dom sequential updating of the sites. The obstacle is 
introduced by modifying the hopping probability to rp 
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at one specific bond along the chain. < r < 1 repre- 
sents a slow bond and r > 1 a fast bond. We choose open 
boundary conditions with the special bond in the middle 
of the chain. 

Mean- field theory predicts an infinitely long traffic jam 
for all r < 1 and only a logarithmic depletion density pro- 
file near a fast bond [21] . The literature is confused about 
what happens in reality. The notion of an rc < 1 queu- 
ing transition seems to have been implicitly presumed 
in the ASEP literature for over a decade. Kandel and 
Mukamel [22] presented early numerical data (for a differ- 
ent but related growth model involving parallel updating 
and polynuclear growth) suggesting a critical point which 
would correspond to « Tc — 0.7 in our model. Above the 
queuing transition, they reported evidence of continu- 
ously varying exponents in the density profiles. Some of 
these aspects have been confirmed by more recent stud- 
ies [23,24], but to the best of our knowledge, the exis- 
tence of an Tc < 1 has not been resolved unambiguously. 
For example, in the ASEP studies with a slow bond and 
periodic boundary conditions, Janowsky and Lebowitz 
presented their phase diagrams as if rc = 1 (the mean- 
field location) [12]. In retrospect Tc — 0.8 seems to us 
consistent with, e.g., their series expansions [24]. Their 
main focus was however elsewhere, with the fluctuations 
in the position of the shock front of the queue, at the far 
end from the slow bond, and not with the precise value 
of rc. 

The exact form of the stationary state in the ASEP can 
often be obtained from the so-called matrix and Bethe 
ansatz method [2,13,15-17], but these analytic techniques 
typically work only for very specific boundary conditions 
and/or update rules. Schiitz [13,25], for example, found 
Tc = 1 for periodic boundary conditions and parallel up- 
dating. That is consistent with our results, because par- 
allel updating creates intrinsically weaker stochastic noise 
than random sequential updating [26]. 

The directed polymer community was focused on the 
slow-bond issue in the mid 1990s [27-32]. The driving 
force behind these studies was the realization of such di- 
rected polymers in terms of flux tubes in type-II dirty su- 
perconductors. The (l+l)-dimensional ASEP is equiva- 
lent to (l+l)-dimensional KPZ-type growth, and the lat- 
ter to a directed polymer in two dimensions subject to a 
random potential. In these equivalences, the slope of the 
KPZ interface is the deviation of the local density from 
a half filling in the ASEP, dh/dx = 1 - 2p{x), and the 
mapping of the KPZ equation to the directed polymer 
problem involves the celebrated Hopf-Cole transforma- 
tion, W = exp[(A/2i/)ft,] (with A and v the KPZ couphng 
constants; for a review see Ref. [32]). The slow bond 
transforms into a columnar defect with a short-ranged 
attractive interaction, and the queuing issue translates 
into whether the polymer becomes localized to it imme- 
diately or only beyond a critical defect strength. The 
latter is true above a critical dimension Dc- Power count- 



ing in the KPZ equation and associated field-theoretical 
renormalization studies suggest that Dc = 1, i.e., our 
ASEP model is at the critical dimension. In such cases 
one expects that Tc = 1, likely accompanied by essen- 
tial singularities [27-32]. Our results presented below 
seem to contradict these field-theoretical studies, but ac- 
tually only do so in a limited sense. We find a more 
complex structure. The queued ASEP phase represents 
the strongly localized state. It exists only beyond a criti- 
cal defect strength rc < 1. The power-law shaped profile 
that remains for weaker slow bonds, represents a form 
of weak localization. Earlier numerical studies in the di- 
rected polymer representation confirmed localization in 
£) = I^c = 1 for all r < 1 [27,28] but were likely insensi- 
tive to this distinction. 

Faced with the realization of this process in terms of 
slow combustion of paper, our first goal is to settle the 
location of rc for the sequential update rule, numerically 
as accurately as possible. We demonstrate here the pres- 
ence of a dynamic phase transition at rc — 0.80(2). The 
queue remains infinite in length all the way up to rc- Its 
density = ^(1 + Ab) decreases as Ab ^ ]rc — r]'' with 
/? ~ 1.46(4). This is presented in Sec. Ill after a detailed 
discussion of our choice of boundary conditions in Sec. II. 

Having settled the existence of the critical point, we 
turn, in Sec. IV, our attention to the density profile near 
the slow bond, p{x) = i[l -|- A(x)], with x the distance 
from the special bond. It follows always a power law, 
A(x) ~ Aft -I- Ax^'^ . Below the transition, r < Tc, the 
density profile has a power-law tail with exponent = ^. 
Above the transition, rc < r < 1 and A;, = 0, the den- 
sity profile has a power-law shape with exponent = ^■ 
Notice that the remaining power-law queue above rc is 
still infinite in magnitude, since jA{x)dx ~ Ns^^ di- 
verges. In the fast-bond scenario, the queue has always a 
power-law profile (i.e., A(, = 0) for all values r > 1, with 
exponent u ~ 0.63(3). 

These power-law density profiles are very intriguing, 
in particular, the fast-bond one. They are different from 
the density profiles near reservoirs, e.g., those in the exact 
solution of Dcrrida et al. [16] for the open boundary con- 
ditions with two reservoirs. Those have exponential tails 
in the reservoir-dominated phases, and power-law tails 
with exponent i in the bulk-dominated maximal-current 
phase. It is well known how to explain these reservoir- 
related profiles with the help of simple scaling arguments 
involving the dynamic exponent z = | , the roughness ex- 
ponent X = ^, and the absence or presence of a nonzero 
group velocity for fluctuations (see, e.g., Ref. [11]). 

In Sec. VI, we generalize these heuristic arguments to 
the density profiles near the slow bond. This reproduces 
the observed slow-bond values = | and u = ^. A 
similar line of reasoning for the fast-bond profile yields 
u = 1/z = ^, but is on shaky grounds, in particular, 
in the light of the fact that our numerical results give 
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systematically a somewhat smaller value. 

Another important result of this study, presented in 
Sec. V, is that the passage of particles through the slow 
bond sets itself up in the stationary state as an uncor- 
related process both below and above Vc- Fluctuations 
travel away from it from both sides. No information 
passes through the slow bond. The fast bond, on the 
other hand, acts very much like a normal site, and fluc- 
tuations flow through it. In Sec. VII, we summarize our 
results. 

II. BOUNDARY CONDITIONS 

The choice of boundary conditions is important for the 

accuracy of our numerical analysis. The special bond 
creates a power-law shaped density profile. Therefore we 
need to fully control other sources for density profiles 
and minimize interference. Wc have chosen open bound- 
ary conditions with the special bond halfway along the 
road, and particle reservoirs on either side, at a; = 1 and 
X = Ns- Particles can only hop to the right, x ^ x + 1. 
The hopping probability is equal to p for all sites, except 
for three special sites. The probability to hop through 
the special bond is equal to p' = rp. The probability 
to enter (leave) the road from (into) the reservoir at site 
a; = 1 (a; = A^s) is equal to ap {(3p). It is advantageous to 
set p = 1, if possible, to maximize the speed of the Monte 
Carlo (MC) simulations, but in our case that would ex- 
clude us from addressing the fast-bond scenario r > 1. 
We set P = \ throughout this study. There are several 
alternatives that can speed up the simulations, but we 
did not feel the need to explore them in this study. For 
an example, one might set p = \ everywhere along the 
chain except at the special bond by increasing the update 
probability of that bond. 

Our choice to employ open boundary conditions might 
appear surprising. They often introduce edge effects (sur- 
face critical phenomena) that are typically more difficult 
to interpret and control than those for periodic bound- 
ary conditions. In the ASEP, however, periodic boundary 
conditions introduce a shock wave in the density profile 
at halfway around the chain, opposite to the slow bond. 
Janowsky and Lcbowitz [12] studied the fluctuations in 
the position of this shock wave and found for it to fluc- 
tuate critically, implying the absence of a characteristic 
length scale. We like to decouple the slow bond from 
those fluctuations and do so by using open boundary con- 
ditions. The trade-off are density profiles near the edges 
of the road induced by the particle reservoirs. But these 
are under full control with the help of the exact solution 
of the r = 1 ASEP with open boimdary conditions [16]. 
The density profiles near the edges have only exponential 
tails, provided that we choose suitable values for a and 
(3. At r = 1, for a = /3 = ^, the density profile is com- 
pletely flat and featureless; px = \ for all x [16]. Thus, 



we select a = (3 = ^ throughout this study [33] 

Assume that the slow bond creates an inflnite queue. 
In the bulk of that queue, far from both the slow bond 
and the road edge, the stationary state is uncorrelated, 
because locally it is indistinguishable from a setup with 
periodic boundary conditions without the slow bond. 
Moreover, from the perspective of the sites near the a; = 1 
edge, this situation is indistinguishable from a setup with 
r = 1 (no slow bond) where an exit probability /3 at the 
opposite site of the road could be responsible for this en- 
hanced bulk density A^. From the exact solution of that 
setup, we know that the density proflle near the entry 
edge is exponential, p{x) ~ i(l -I- A;,) -I- i3cxp(— .t/^) 
with a finite correlation length ^ ~ A^. Our simulations 
confirm this [34]. 

Other important features of the density profile arc pre- 
determined as well. For all a = (3, the density profile has 
particle-hole symmetry with respect to the special bond, 
p{x) = 1 — p{Ns + 1 — a;). Define A(a;) as the deviation 
of the density from i, p{x) = ^[1 -|- A(a;)], such that 
A(x) = —^{Ng -|- 1 — a;). Then, Ab is the order parame- 
ter of our model and represents the spontaneous faceting 
angle of the slow-combustion interface profile in the KPZ 
interpretation. 

In the steady-state limit, the current along the chain 
must be uniform. From the fact that the bulk station- 
ary state is uncorrelated, it follows immediately that its 
value through such a bulk bond is equal to 

J = p(n,(l-n,+i)) = |(l-A2). (1) 

This means that we have the option to determine A(, by 
measuring J. The current from the reservoir to the first 
site 

J = ap((l-ni)) = ^(l-Ai) (2) 

must be equal to the bulk current in the stationary state. 
This yields, for a = ^, that the density at site a; = 1 is 
equal to Ai = A^. 

The density profiles in front and beyond the special 
bond obey particle- hole symmetry, (til) = 1 — {nn), with 
XL and xn the sites immediately in front and beyond the 
special bond. Moreover, we will demonstrate in Sec. V 
that the ratio IZ = {fiLfiR) / {nn) is very close to 7^ = 5 
for all values of r. The current through the special bond 

J = rp{hL{l - ur)) = rp [(1 + 7^)(ni) - U] (3) 

must again be equal to the bulk current, thus anchor- 
ing the value of the density immediately in front of the 
special bond to the order parameter of our process as 
Al = |[(1 - Ag)/r - 1] if 7e is exactly equal to \. 

Figure 1 summarizes the above discussion. The cur- 
rent, the densities at the first and last sites, the char- 
acteristic exponential length scale of the density proflle 
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near the reservoir edge (in the faceted phase), and the 
bulk density are all linked to each other. This leaves only 
the spontaneous creation of a nonzero A;, and the density 
profile near the special bond as independent issues. 

III. QUEUING PHASE TRANSITION 

The first issue at hand is to settle by numerical means 
the burning question whether there exists a queuing tran- 
sition at an Tc < 1. We perform MC simulations for sys- 
tem sizes up to Nf, = 4096 and analyze them by finite 
size scaling (FSS) techniques. The order parameter of 
the queuing transition is the offset of the density in the 
bulk, Afc (far from both the special bond and the edge). 
There are several ways to measure this. 

One can directly measure the density at a site such 
as a; = jN^ and perform a FSS analysis to determine 
the asymptotic value. This works, but neither x — jNg 
nor any other fixed site is optimal for such an analysis, 
because the density profile has a power-law tail at the 
special bond side and only an exponential one near the 
reservoir. Instead, we show in Fig. 2 the values for 
from a power-law density profile fit 

A{x) ~ Ab + Ai-" (4) 

(with X = XR — X the distance from the special bond) at 
our maximum system size Ns = 4096. 

As pointed out in the preceding section, Af, is directly 
linked to various other quantities, such as the average 
current J, the density at the first site near the edge Ai, 
and the characteristic length ^ of the exponential tail 
in the density profile near the edge. We measure these 
quantities as well, and translate the first two into their 
predictions for Aft. The results, also shown in Fig. 2, 
arc almost indistinguishable from those of the power-law 
density profile fits. This confirms our analysis of the pre- 
ceding section. 

Figure 2 suggests very strongly the existence of a crit- 
ical point at about Tc — 0.8. However, this is a common 
optical illusion, which vanishes upon zooming-in to this 
point, as in Fig. 3. Assume that the order parameter 
obeys the conventional FSS scaling form 

AbiN,,e) = 6---Afe(6-i7V„ 6^e), (5) 

where e = Tc — r. We test how well our numerical data 
obey this scaling relation and what the best values of Vc, 
xa, and P = xa/v are. The order parameter should 
scale as a function of e as A;, ~ . In Fig. 4, we 
show a log-log plot of A^ versus e for various choices 
of at iV = 4096. The best straight line is obtained 
for rc = 0.80(2) with slope (3 = 1.46(4). At the same 
choice for rc, the order parameter also scales perfectly 
as a power law A;, ~ iV~^^ with a critical dimension 



XA = 0.370(5). One should always be on guard for cor- 
rections to scaling. For that reason we plot in Fig. 5 the 
scaling function <S defined as 

At,{N,,e)=N-^-S{Nye) (6) 

for Tc = 0.80, XA = 0.370, and p = xa/v = 1-46. The 
data collapses very well, implying only minor corrections 
to scaling. 

An alternative scaling form to consider is an exponen- 
tial esscntial-singularity-type infinite-order transition, in 
particular, with Vc = 1, as suggested by the directed 
polymer renormalization studies [27-32] . We tried these 
forms, shown in Fig. 6. They fit our MC data poorly. 

IV. DENSITY PROFILES 

The density profiles near the special bond have a 

power-law shape for all values of r. Figure 7 shows our 
numerical results for the exponent v and the amplitude A 
as defined in Eq. (4). We performed also two-parameter 
fits after determining Af, independently from the numer- 
ical values for the current and the density at site x = 1, 
using the inter-relations outlined in Sec. II. These results 
are identical within the MC noise. 

The jumps in at Vc and r — 1 are very pronounced in 
Fig. 7. It seems safe to conclude, surely as a starting as- 
sumption for the discussion in the following two sections, 
that for slow bonds the exponent takes the value v — ^ 
in the r < Tc macroscopic queued phase and = | in the 
Tc < r < 1 power-law queued phase. To the best of our 
knowledge only one earlier study, the one by Slanina and 
Kotrla [23], observed this type of power law, but they 
suggested a value different from v = Wc will present 
convincing heuristic analytic derivations for our values in 
Sec.VI. 

The power-law queue for fast bonds, r > 1, is quite 
intriguing. This is where Kandel and Mukamel [22] 
sighted a possible continuously varying v. We interpret 
our data as strong evidence for a nonvarying constant 
value V ~ 0.63(3). The drop in the estimates in Fig. 7 
near r = 1 resembles conventional (multicritical-type) 
crossover scaling, but we cannot verify this explicitly, 
because this power-law decays much faster than in both 
slow-bond phases, and, e.g., at r = 1.1, the amplitude 
sinks underneath our MC noise level already at about 
x ~ 60. We exclude a: < 20 from our fits to avoid (short- 
distancc-type) corrections to scaling. 

Obviously we would like to "talk" this fast-bond 
power-law profile towards p = ^, since that number oc- 
curs naturally in ID KPZ-type processes. However, our 
heuristic argument for = 1/z = | is not very strong, 
see Sec.VI, and the fits to the MC data in Fig. 7 remain 
consistently below that value. 

For this reason we also studied the following r ^ oo 
like setup. Consider a normal chain without any special 
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bond but with the site in the middle aUowed to be doubly 
occupied, nMj2 = 0,1,2. In this setup we can increase 
the hopping probability p to p = 1 and thus speed up 
MC simulations {a and (i arc again set equal to |). The 
log-log plot of the density profile, shown in Fig. 8, is quite 
straight. Still, the slope suggests a somewhat smaller ex- 
ponent, V = 0.64(2) (using 20 < .t < 500 as the fitting 
range), i/ = | is still a possibility, but seems to require a 
significant subdominant correction to scaling power-law 
term. A more detailed analysis becomes meaningful only 
when the noise level is brought down by at least one more 
order of magnitude from our current 5p/ p ~ 0.001 level, 
which requires vastly longer MC runs. 

V. UNCORRELATED PASSAGE 

Our numerical observation that the density profile near 
the special bond follows always a power law, 5p ^ x~'^, 
in all three phases, is far from obvious. The actual values 
for v are even more intriguing. In this and the follow- 
ing sections, wc present intuitive heuristic explanations 
for the slow-bond values, and also address the fast-bond 
case. An important ingredient in this is that the passage 
through the slow bond is an imcorrclatcd random process 
in both the macroscopic queued phase at r < rc and the 
power-law queued phase at rg < r < 1. 

In the macroscopic queued phase, the absence of pas- 
sage correlations is easily understood. Fluctuations 
travel away from the slow bond, both in front and be- 
yond it. The group velocity of fluctuations Vg points 
away from the slow bond in both directions. Vg = 6 J/ dp 
represents the local response of the current to a density 
fluctuation. The stationary state is uncorrolatcd inside 
the bulk, such that the current is equal to J = pp(,(l — p(,) 
and Vg = p{l — 2pb) = —pAb- Fluctuation-type wave 
packets travel with this velocity along the road, while 
they broaden spatially as ^ ~ t^/^, with the ID KPZ 
dynamic exponent -2 = §• In the KPZ growth context, 
Ab represents the average slope of the growing surface, 
and the traveling wave packet reflects that the interface 
moves perpendicular to the local surface orientation. 

The precise form of Vg for spatially varying densities 
p{x) is more complex, but for slowly varying ones, like 
here, we can assume Vg is well represented by 

Vg{x)=p[l-2p{x)] = -pA{x). (7) 

A(x) is positive in the macroscopically queued phase, 
such that the center of mass of a fluctuation packet moves 
away from the slow bond linearly in time, xqu ~ t- Dur- 
ing this process, it spreads over a width ^ ~ t^^^. Fluc- 
tuations detach from the slow bond, because the center 
of mass of the packet propagates faster than its broaden- 
ing front. Therefore the density fluctuations at the slow 
bond are uncorrelated in time. No memory remains at 



the slow bond of anything happening there before. No 
information passes through the slow bond. 

Most of this remains true in the power-law queued 
phase at Tc < r < 1. Now the group velocity van- 
ishes in the bulk, but remains nonzero near the slow 
bond, because of the power-law shaped density pro- 
file of the queue. The center of mass of a fluctuation 
packet still moves away from the slow bond, but only as 
xcM ~ t^^^^~^''\ see Eq. (7). During this, it spreads again 
over a width ^ ~ t^^^. Therefore, for all < 2 — 1 = ^ the 
packet detaches from the slow bond, p = z — lis the crit- 
ical value. Numerically we find i' — ^, see Fig. 7(a). So 
the density profile near the slow bond organizes itself in 
a form where the passage fluctuations through the slow 
bond are uncorrelated in time and density fluctuations 
originating on the road do not affect it. 

It is useful to test this explicitly by numerical simula- 
tions, in particular, in the power-law queued phase. Time 
correlators, such as the current-current autocorrelation 
function, are the preferred tools for this, but unfortu- 
nately they do not yield much useful information. The 
current-current correlator drops in magnitude by two or- 
ders within 10 MC time steps, not only near the slow 
bond, but everywhere along the road as well. This re- 
flects that, in KPZ growth, 5 J ~ N~"' scales with a large 
cr ~ 2. 

As a second best choice, we focus instead on spatial 
correlations between the densities across the slow bond. 
Consider the ratio 

7^ = %^, (8) 

with fiL and hji the density operators at the sites imme- 
diately in front and beyond the special bond. For refer- 
ence, the same type of ratio for two nearest-neighbor sites 
anywhere along the road, in the bulk, far from edges and 
slow bonds, is equal to i, because the bulk stationary 
state is uncorrelated, with (nin^+i) = (hi){hi^i). Near 
the edges, however, and in particular, inside power-law 
profiles, the neighbors are correlated, and the ratio moves 
away from ^. 

Figure 9 shows that TZ, as defined in Eq. (8), is al- 
most equal to ^ for all values of r. The deviations are 
only of order 3%. This is consistent with the picture that 
fluctuations travel away from the slow bond from both 
sides. 

To quantify this in more detail, we consider the follow- 
ing mean-field-type approach, in which the road in front 
and beyond the slow bond are treated as reservoirs (de- 
void of fiuctuations as far as the slow bond is concerned). 
We solve thus the following two-site problem, with only 
sites xl and xr on either side of the slow bond. Particles 
hop onto site Xl with an effective probability a^ftp from 
the road in front of it, treating the road as a reservoir; 
then move through the slow bond with probability rp; 
and finally hop away from xr onto the road beyond the 
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slow bond with probability acsP, treating that as a reser- 
voir as well. Finally, we tune a^sp to the value where the 
current takes the same value as in the true system. This 
approximation yields rather trivially TZ = ^ for all r and 
all aeff- The dashed line in Fig. 9 shows TZ for the next 
level of mean-field theory, with four sites instead of two, 
taking into account local correlations. These suffice to 
reproduce already most of the small deviations we ob- 
serve in the true 7?. as a function of r, and support the 
uncorrclatcd passage nature of the process. 

The ratio remains equal to 7^ ~ ^ for fast bonds, 
and actually even better than for r < 1. How do fluc- 
tuations travel there? The group velocity changes its 
sign (moving direction), because the power-law queue, 
A{x) ~ Ax~" , turns into a depletion zone with negative 
amplitude A. Fluctuations travel towards the fast bond 
from both sides. It might seem therefore that this pas- 
sage process must be highly correlated. However, fluc- 
tuations originating from all over the road bombard the 
fast bond from both sides, and average each other out. 
Consider a fluctuation created at a distance x from the 
fast bond. The center of mass of this fluctuation moves 
towards it, and arrives after a time of flight t ~ x^^^ . 
During this time, it has broadened over a width ^ ~ f^/^. 
Ignoring the center of mass movement, its leading edge 
would arrive at the fast bond after a time f ~ a;^. For 
u > z — 1 = \ the leading edge arrives well before the 
center of mass, and the latter can be neglected. Again, 
V = z — 1 = ^,\s the critical value. For fast bonds we flnd 
numerically v ~ 0.63, see Fig. 7(a). This explains why 
TZ = ^. The density proflle organizes itself again into a 
form where the passage correlations remain simple. The 
fast bond acts very much like an ordinary bulk site, and 
fluctuations flow through it. 



VI. A DERIVATION OF THE DENSITY 
PROFILES 

In the preceding section we found that the density pro- 
files organize themselves into a form such that the pas- 
sage through the special bond is an uncorrelated pro- 
cess. Here we give heuristic arguments for the actual 
values of the exponents: v = j at r < and u = ^ at 
Tc < r < 1 for slow bonds. We also explain why for fast 
bonds z/ ~ 1/z. 

First, consider the slow bond = power-law density 
profile in the macroscopic queued phase at r < Tc. The 
total number of excess particles in this queue diverges 

1 /2 

with system size Ng as SN ^ Ns . This has a familiar 
ring to it. In the Derrida et al. [16] type open system 
setup with reservoirs on both sides and no special bonds, 
the fiuctuations in the total number of particles on the 

1/2 

road scale as SN ^ Ng' . In that setup, this property 
does not translate into power-law-type density profiles. 



except when the road is half filled, ph = ^. The density 
profiles are exponential or featureless in the two Pb ^ \ 
phases where either reservoir controls the bulk density. 

The parking garage process of Ref. [11] is closer to 
queuing dynamics. In that study the two reservoirs were 
merged into one, such that the road forms a loop start- 
ing and ending in the same parking garage. The total 
number of cars in the system is then conserved, leading 
to dynamic phase transitions between condensate-type 
stationary states where the garage is macroscopically oc- 
cupied and a normal phase where it is not. That process 
has two parameters, the total number of cars in the sys- 
tem and a modified hopping probability ap to jump from 
the garage onto the first site of the road. 

The fiuctuations in total number of particles on the 

1/2 

road is again equal to 5N ^ Ng . The explanation 

of this goes as follows [11] for the normal phase, the 
nonmaximal-current condensate phase, and also at the 
transition point between them. The group velocity of 
fluctuations Vg is nonzero. This means that the depar- 
ture of cars from the garage is an uncorrelated process. 
Fluctuations detach from the garage because they travel 
away faster (linear in time) than they arc spreading back- 
ward (as ^ ~ f-^/^ with z = |). Moreover, after a time of 
flight tflight = Ng/Vg they move around the loop, return 
to the garage, and arc completely erased. So we deal with 
^flight random uncorrelated deposition events. The fluc- 
tuations in the number of cars on the road therefore scale 

1/2 

as tf/igjji^. In the condensate phase, these fiuctuations do 
not lead to an offset in the average density of parked cars, 
because the garage is macroscopically occupied, and pos- 
itive and negative fiuctuations cancel out against each 
other. But at the transition point, the bottom of the 
garage becomes visible. This limits the negative density 
fiuctuations, and therefore introduces a bias towards in- 
creased occupation, such that the number of parked cars 
is enhanced and scales as 5Np ~ nI^"^ . 

The same type of reasoning applies to the slow-bond 
setup. The passage through the slow bond is a stochas- 
tic uncorrelated event (as demonstrated in the preceding 
section), similar to departures from the garage mentioned 
above. Again, all memory is erased after tflight ~ Ng^ the 
time that takes for a fluctuation to travel from the slow 
bond to the reservoir. The fiuctuations in the number of 

1/2 

cars passing through therefore scale as ifj-ght- These fluc- 
tuations are biased again, because the sites immediately 
in front and beyond the slow bond are not reservoirs. 
Excess particles waiting to pass are spread out, and not 
available for immediate passage. The passage process is 
biased, because slow bonds process particles slower than 
normal bonds, while the passing probability of vacan- 
cies does not depend on the value of r. (Our process has 
particle-hole symmetry but only in conjunction with left- 
right mirroring with respect to the special bond.) The 
total number of excess cars near the slow bond waiting to 
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1/2 

pass scales therefore as SNp ^ Ns . These extra parti- 
cles must be accommodated over a stretch of road x < xl 
behind the slow bond. We can imagine two ways to real- 
ize this: an exponential density profile with a correlation 

1 /2 

length diverging as ^ ~ A^g ' or a power- law profile with 
u = ^ as we actually observe. The power law is indeed 
more likely given the intrinsic critical nature of ASEP. 

Next, let us generalize this argument to the p = ^ 
power law above the queuing transition, at rc < r < 1. 
The time of flight of a fluctuation to travel from the slow 
bond all the way back to site x = 1, scales now only as, 
iflight ~ since Vg = —pA{x) ~ x~''. Assume that 

< ^, in which case the fluctuations still detach from 
the slow bond and the passage through the slow bond re- 
mains uncorrelated. The process at the slow bond is still 
biased toward low density fluctuations. FSS corrections 
to the total number of cars in the queue is proportional 
to fflight uncorrelated events: 

6N ~ ^ N(r^y\ (9) 

This queue heaps up behind the slow bond, and again 
arranges itself in the form of a power-law shaped den- 
sity profile, 6p x~''. Self-consistency implies that 
(v + l)/2 = —V + 1 ^ V = ^, m accordance with the 
observed value. 

The V ~ 0.63(3) power law for the fast bond is 
more challenging. This is a fundamentally different phe- 
nomenon. Again fluctuations travel across the system, 
but now run towards the fast bond instead of away from 
it. Actually, as shown already in the preceding section, 
for J/ > z — 1 = i the time of flight of the center of mass 
of a fluctuation ifiight ~ Ns~^^ is longer than the time 
it takes that same fluctuation to spread over the entire 
system t ^ N^. In the r < 1 phases, we are allowed to 
ignore for this reason the spreading of the fluctuations, 
and only consider their center of mass motion (the time 
of flight). At r > 1 this is reversed. In the r < 1 phases, 
the exponent v was insensitive to the actual value of the 
dynamic exponent z of the dynamic process. In the r > 1 
phase, it must depend on z. 

A ~ I power-law tail is very rare. It does not 
appear, e.g., anywhere in the Derrida et al. [16] type 
two-reservoirs setup. Interestingly, however, this den- 
sity profile appeared already in the parking garage ASEP 
study [11]; at the second- type condensation transition, 
from the "normal" to the "maximum current" phase. 
The characteristic feature was that, at the transition 
point, the garage started to transmit information (seized 
to act as a reservoir), and that at that point the bulk 
group velocity was zero. The similarities with the fast 
bonds are striking. We are clearly looking at the same 
type of phenomenon. The slow bond does not trans- 
mit information, while the fast bond acts very much like 
a normal site and fluctuations move (i.e., they spread) 
through it (see also the preceding section). 



What might the true value of ~ 0.63(3) be? An ob- 
vious guess is that u = 1/z, but how to explain this? One 
of the crucial aspects must be again that the processing 
of fluctuations is biased at the fast bond, leading to a 
depletion queue with a total deficit of 

rNs/2 

Jo 

particles. The sign of this in now negative, because at the 

fast bond the particles are processed faster than at the 
normal bonds, while the passing rate of vacancies does 
not depend on r. 

Density fluctuations are created everywhere along the 
road, all the time, and with a common characteristic am- 
plitude. Each spreads in time over a region ^ ~ t^/^. A 
fluctuation created at a distance x from the fast bond 
arrives there after a time t ~ a;^ and with a reduced 
amplitude (from spreading) of order Ax~^^''. The asym- 
metry in processing high and low density fluctuations 
gives rise then to a density deficit of order Ax"^^^ from 
fluctuation originating at distance x. Next, adopting 
rather frivolously superposition principle concepts, one 
would guess that the total density deficit in front of the 
fast bond scales as SN ~ j^"^'^ dx Ax~^/^ ~ Ns~^^^ in 
agreement with what we observed numerically. 

Although the last argument is reasonably appealing, it 
is certainly not convincing. A more robust explanation 
is needed. Moreover, the numerical value i/ ~ 0.63(3) 
is sufficiently lower to cast serious doubts that = | is 
correct. On the other hand, the above argument serves 
as a proper order of magnitude estimate. 

VII. SUMMARY 

In this study, we reconfirmed numerically, and beyond 
doubt the presence of a queuing phase transition in the 
ASEP with a slow bond of strength rc ~ 0.80(2). We 
established the two scaling exponents of this transition. 
The order parameter, the excess density in the queue, 
vanishes as Ah ~ \rc — '"['^ with /3 = 1.46(4). At the tran- 
sition point, the number of particles in the queue scales 
with system size as ^ N^^ with xa = 0.370(5). 

From a more general perspective, the transition illus- 
trates that weak obstructions do not give rise to macro- 
scopic traffic jams (queues with lengths that scale linearly 
with the system size). The stochastic fluctuations over- 
whelm the slow bond above rc- 

A second result of our study is that above rc a power- 
law shaped queue (traffic jam) remains, 5p ~ Ax"", with 
X the distance from the obstruction and opposite signs 
for A in front and beyond the obstruction. The exponent 
is equal to v = |. This value is most likely universal, be- 
cause our heuristic derivation for the general case does 
not involve specific details of the dynamics. In partic- 
ular, it does not involve the KPZ value of the dynamic 
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exponent z, except for the requirement that fluctuations 
travel faster away from the slow bond than they spread. 
The argument applies for any process with z > l + v = ^. 

From the directed polymer perspective our results are 
unexpected. The slow-bond queuing transition repre- 
sents a crossover from strong to a weaker (but still strong) 
form of localization, because the >^ = \ power-law density 
tail near the slow bond contains still an infinite number 
of particles; (naively) the polymer distribution behaves 
as {W) ~ exp[(A/2i^)(/i)] - CTcpi-Cx^-"). It will be in- 
teresting to see how this weak localized phase, and the 
above exact self-consistent argument for z/ = |, can be 
integrated and reconciled with the field-theoretical de- 
scriptions of Refs. [27-32]. 

For fast bonds (such as a local widening of the road), a 
macroscopic depletion queue (with a length proportional 
to the road length) never appears. Instead, a power-law 
shaped depletion queue is always present with exponent 
V ~ 0.63(3). It remains yet imclear whether this value is 
equal to v = 1/ z. It will be interesting to study how our 
results extend to other models of driven fiow along one- 
dimensional channels, in particular, to non-KPZ-type dy- 
namics. 
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FIG. 1. Schematic density profiles are shown: for the slow 
bond (SB) (a) r < (queued phase) and (b) rc < r < 1 
(nonqueued SB phase) and for the fast bond (FB) (c) r > 1 
(nonqueued FB phase). 
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FIG. 2. The order parameter Ab vs the strength r of the 
special bond at Ns = 4096, determined from three different 
datasets: the average current J (squares), the density Ai at 
the first site neax the reservoir edge (crosses), and three par 
rameter power-law fits to the density profiles near the special 
bond (circles). 
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FIG. 3. The order parameter Ai, vs the strength r of the 
slow bond in the vicinity of the critical point rc, as obtained 
from the average current dataset. 
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FIG. 4. Determination of the critical point and critical ex- 
ponents, (a) Double logarithmic At ~ |e|'' type plots of the 
order parameter with e = — r at Ns = 4096 for various 
choices of rc. The best straight line is found at rc = 0.80(2) 
and with slope /3 = 1.46(4). (b) Double logarithmic plots of 
Ai, ~ N^^^ as a function of system size Ns at rc = 0.80. 
The slope (dashed line) yields xa = 0.370(5). For clarity we 
show only the data for A^ obtained from J (squares) and Ai 
(crosses) . 
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FIG. 5. Data collapse by Eq. (6), i.e., the scaling function 
of the order parameter using the values rc = 0.80, xa = 0.370, 
and /3 = 1.46 as found in Fig. 4 
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FIG. 6. The same data as in Fig. 2 fitted to a scaling 
form of the type A;, (A;, = exp[— o(l — r)**]) represent a 
so-called essential-singularity characteristic for a possible in- 
finite-order-type transition with rc = 1. The curves fail to 
straighten out, indicating this is a poor fit. 
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FIG. 7. The exponent ;/ (a) and the amplitude A (b) of 
the power- law shaped density profile, defined in Eq. (4), for 
various r at system sizes Ng = 2048 (crosses) and Ng = 4096 
(circles) . 
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FIG. 8. The density profiles in the r ^ oo model at 
Ns = 4095, implemented as a normal chain with uniform 
hopping probability p = 1, but one special double occupancy 
site in the middle. The dashed line, with slope v = 0.64, 
serves as guide to the eye; a slope v = 2/3 seems too steep. 
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FIG. 9. The ratio TZ defined in Eq. (8). TZ remains close 
to the uncorrelated passage value i for all r. The small de- 
viations, of order 3%, do not scale with system size, and axe 
mostly described already by a four-sites-type mean-field ap- 
proximation, the dashed line. 
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